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Many real-world networks describe systems in which interactions decay with the distance between 
nodes. Examples include systems constrained in real space such as transportation and communi- 
cation networks, as well as systems constrained in abstract spaces such as multivariate biological 
or economic datasets and models of social networks. These networks often display network motifs: 
subgraphs that recur in the network much more often than in randomized networks. To understand 
the origin of the network motifs in these networks, it is important to study the subgraphs and 
network motifs that arise solely from geometric constraints. To address this, we analyze geometric 
network models, in which nodes are arranged on a lattice and edges are formed with a probability 
that decays with the distance between nodes. We present analytical solutions for the numbers of 
all 3 and 4-node subgraphs, in both directed and non-directed geometric networks. We also ana- 
lyze geometric networks with arbitrary degree sequences, and models with a field that biases for 
directed edges in one direction. Scaling rules for scaling of subgraph numbers with system size, 
lattice dimension and interaction range are given. Several invariant measures are found, such as 
the ratio of feedback and feed-forward loops, which do not depend on system size, dimension or 
connectivity function. We find that network motifs in many real-world networks, including social 
networks and neuronal networks, are not captured solely by these geometric models. This is in 
line with recent evidence that biological network motifs were selected as basic circuit elements with 
defined information-processing functions. 

PACS numbers: 05, 89.75 



I. INTRODUCTION 

Many systems in nature can be represented as com- 
plex networks 0, 0, 0, H- Often, natural and engi- 
neered networks have an underlying geometric arrange- 
ment. In such networks, nodes are embedded in a ge- 
ometric space, and edges tend to link nodes that are 
close neighbors. Examples include the physical layout of 
the internet JEJE IE ff > transportation networks jjj 0] , 
power grids as well as networks of wiring between 

neurons 0, Il2lll3l IkH ] or cortical areas 0, 0] . 

A geometric constraint does not have to be of spatial 
origin. In models of social 17] or world-wide-web net- 
works 0, 0], nodes may be assigned attributes (e.g. 
language and field of interest of web pages, occupation 
and residence in social networks) and links may be corre- 
lated with closeness in this attribute-space [2(J, |2j| . Geo- 
metric networks can also arise from analysis of multivari- 
ate datasets. For example, networks have been proposed 
to describe gene expression (2^ or economic datasets [23| , 
where distance between nodes corresponds to high corre- 
lation coefficient in the dataset. These networks can be 
embedded in a high dimensional Euclidean vector space. 

Many of these diverse networks have been recently 
found to display network motifs 0, 0, 113 : a charac- 
teristic set of recurring subgraphs which occur much more 
often than in randomized networks with the same degree 
sequence. Different networks usually display different 
motifs, and motifs can be used to characterize families 
of networks. In biological networks, it has been demon- 
strated that each of the motifs can carry out a key infor- 



mation processing function [H EHH |H HJHJ H| • 

It is of interest to study the origin of network mo- 
tifs in each real- world network. In particular, it is of 
interest to compare the network local structure to that 
of model networks which have a similar geometric con- 
straint. For example, the researchers who mapped the 
synaptic wiring of C. elegans speculated that :"The abun- 
dance of triangular connections in the nervous system of 
C. elegans may thus simply be a consequence of the high 
levels of connectivity that are present within neighbour- 
hoods" 0]. An analysis of the abundance of subgraphs 
in purely geometric networks, can help discern whether 
the motifs arise based on simple geometry, or whether 
they arise due to additional optimization or design of the 
network [33ll3^| . 

To address this, we study geometric network mod- 
els, in which nodes are arranged on a lattice, and edges 
are placed randomly between nodes with a probabil- 
ity F(r) that decays with the distance between nodes. 
Several features of related models were previously stud- 
ied P, |M E^LBI ■ These features include degree 
distributions |35ll3q|. diameters [37j and clustering coef- 
ficients . Here we focus on the subgraph content of 
these networks. We consider non-directed and directed 
networks, as well as cases where directed edges are biased 
in a particular spatial direction. We present an analyti- 
cal solution for the numbers of small subgraphs and the 
scaling of all types of subgraphs with system size and lat- 
tice dimensionality. We find invariants that can be used 
to easily compare networks with geometric models. 
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FIG. 1: Non-directed geometric networks and their randomized versions, (a) Id geometric network with N = 40, R = 7, (k) = 2. 
(b) Id randomized version, (c) 2d geometric network with N = 900, R = 1, (k) = 2. (d) 2d randomized version. The 2d network 
shown has non-continuous boundary conditions, for clarity. 



II. RESULTS 

A. Non-directed geometric model 

In the geometric model, N nodes are arranged in a 
(i-dimensional Euclidean lattice with toroidal (continu- 
ous) boundary conditions (Fig 1). Non-directed edges 
are placed at random according to a connectivity density 
function F(x,y) (where F{x,y) = F(r(x,y)), and r is 
the distance between nodes x and y). For each pair of 
nodes, x and y, a random number p is generated, and an 
edge is placed if p < F(x,y). F(r) is a decaying func- 
tion with a range R — J rF(r)dr. We consider the case 
where R is much larger than the lattice spacing [R » 1) 
but where the effective connectivity neighborhood of each 
node is much smaller than the system size R d << N. In 
this case, the mean number of edges per node is 




The degree distribution, the distribution of number of 
edges per node, P(k) is Poissonian, with a mean of 
(k) (assuming that F(r) decays sufficiently rapidly (65|). 
Therefore, a random network ensemble that preserves 
the degree distribution of the geom etric network, is the 
Erdos-Renyi model | 4(1 l4ll |42j, |43j with N nodes where 
edges are placed at random with probability pErd = 
(k)/N. 

We now calculate the mean number of appearances of a 
given subgraph in the geometric model. The probability 
for the subgraph may be expressed in terms of overlap 



integrals of the function F(r). For example, the triangle 
subgraph tends to occur when three nodes, x, y and z are 
sufficiently close, as expressed by the integral 

(G A ) = j j j F(0, y), F(0, z), F(y, z)dydz (2) 

Where without loss of generality x is at the origin. The 
factor 1 /6 is due to the symmetry, where the same trian- 
gle can be counted if y or z serve as the origin, and if y 
and z are interchanged when x is the origin. The sym- 
metry factor can be calculated based on the symmetry of 
each subgraph (one over the number of permutations of 
nodes that lead to an isomorphic subgraph). 

In table I, we present the number of appearances of 
all three and four-node subgraphs in a non-directed ge- 
ometric network. The results in the table apply to the 
case of sparse networks, where (k) << R d (as occurs in 
almost all real world networks). The results are for two 
connectivity functions (Fig 2). The first is a Gaussian 
connectivity function, where 

F g (r) = k(2irR 2 )- d/2 exp(-r 2 /2R 2 ) (3) 

and where 

r 2 = J2(xi-Vi) 2 (4) 

i 

is the L 2 norm, {xi}, {?/.;} denote the d-dimensional co- 
ordinates of nodes x and y. The second connectivity 
function is a hard-cube connectivity function : 

F c (r) = k(2R)- d e(r < R) (5) 
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FIG. 2: Gaussian (dashed) and hard-cube (bold) connectivity 
functions. 



where 



2/0 1 



(6) 



is the L°° norm, and is a step function. Similar over- 
lap integrals (Appendix A) appear in the calculation of 
virial coefficients |57) . and in calculations of percolation 
thresholds lEH. 



Scaling of subgraph numbers with system size 
and dimension 



We present a simple scaling argument for the subgraph 
content of geometric networks. In this picture, the neigh- 
borhood of each node at distances smaller than R is sim- 
ilar to an Erdos-Renyi network. The number of appear- 
ances of a subgraph with n nodes and g edges in an Erdos- 
Renyi network of size Ne an d mean connectivity (k) has 
been shown to scale as |43|, |44| 



(G 



End 



(7) 



In total, there are on the order of N/R d such Erdos- like 
domains in the entire network, each one of size Ne = R d - 
Therefore, the scaling of the number of appearances of 
subgraph G in the geometric network is: 



{Ggeom ) 

which results in : 



(N/R d )N E n ~ 9 {k) & 



(G geom ) ~ NR^-a-^^kf 



(8) 



(9) 



All subgraphs in the geometric network scale with net- 
work size N as G geom ~ N (Table I) . This is in contrast 
to Erdos networks, in which G ~ N n ~ 9 . Therefore, 
in the limit of large system size, all subgraphs in which 
g > n will be network motifs. This includes triangles, 
squares, and aggregates of triangles. 



FIG. 3: Scaling regimes of the geometric model (N = 
5000, d = l,k — 4). At low R/N ratios triangles scale as 
(R/N)^ 1 . At high R/N ratios the network becomes an Erdos 
network and triangles scale as (R/N)°. 



The present scaling argument also provides the regime 
in which finite-size effects begin to be important. Finite 
size effects begin when R d ~ N. When R d > N, the 
entire network essentially behaves as an Erdos netwok, 
and the scaling crosses over to Erods-Renyi scaling (Fig 
3, see also |3q). 

The scaling relations can also be written in terms of the 
network clustering coefficient 0,0iS>0iEEEjE3i which 
is related to the ratio of triangles to V-shaped subgraphs 
(Table I): 



C = 3 
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R d 



(10) 



Inserting this into Eq. (9) we get: 

{G geom (n,g))~ NC ^n+i {k) ^i)- g{d -i) (n) 

In geometric networks, knowing C and (k) is sufficient 
to find the scaling of all subgraphs. This is not the case 
for other types of networks: In general, in non-geometric 
networks C does not determine the number of four-node 
patterns (e.g. diamonds - pattern 6 in Table I) and larger 
motifs. 

C. Directed isotropic geometric model 

We now consider the case where each edge has a direc- 
tion, yielding a directed network. This network is built 
in the same way as the non-directed network, except that 
each pair of nodes x,y is considered twice, and directed 
edges can connect x to y and y to x. The mean number 
of outgoing edges per node is equal to the mean number 
of incoming edges: 



(fcm) = (kout) = / F(r)dr 



(12) 
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TABLE I: Numbers of non-directed three and four-node subgraphs in the geometric model with N nodes, mean connectivity 
(k), range R and dimension d. Pre-factors, /, are for Id and 2d Gaussian connectivity function (fa) and hard-cube connectivity 
function (fhc)- Also shown are the mean number of subgraphs in Erdos networks with mean connectivity (k). Stars represent 
subgraphs that are network motifs in the limit of large system size. 
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TABLE II: Numbers of directed 3-node subgraphs in the geometric model with N nodes, mean connectivity (k), range R and 
dimension d. Subgraphs are grouped according to number of edges. Pre-factors, /, are for Id and 2d Gaussian connectivity 
function (/g) and hard-cube connectivity function (fhc)- Field factors are for the Id Gaussian model (a = 0.27, b — 0.365, c = 
0.73, d = 0.135). Also shown are the mean number of subgraphs in Erdos networks with mean connectivity (k). Stars represent 
subgraphs that are network motifs in the limit of large system size. Note that subgraphs 4, 5, 8 are not network motifs when 
compared to randomized networks that preserve the degree distribution of both single and mutual edges [3, 12a . |44| 
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TABLE III: Subgraph ratios in real world networks and in the geometric model. Networks are (i) social network 1 of 0, A = 67, 
E = 182, (ii) WWW hyperlink network 3 of § (source 13]), N = 47978, E = 235441. (iii) neural synaptic connections in 
C. elegans, N = 280, E = 2170 [lEt (iv) strong neural synaptic connections in C. elegans, N = 280, E = 400 14] (only 
connections with 5 or more synapses). Brackets indicate the absolute count for the subgraph taken as 1 in the ratio. The ratios 
for the isotropic geometric model (p = 1) apply to the poissonian model and for any dimension and connectivity function. 
Ratios 3 and 4 in the isotropic model apply also to geometric networks with arbitrary degrees. The ratios for the fully biased 
anisotropic geometric model (p = 0) apply to the Id geometric model with gaussian connectivity function. In ratios 3, 4 for the 
WWW network a sampling algorithm for subgraph counting was applied (56| with 10 6 samples. 



We find that the same scaling arguments hold. The 
subgraphs fall into classes according to the number of 
nodes and edges n, g. Table II shows the result for three- 
node subgraphs for d = 1 and 2. The 3- node subgraphs 
fall into five classes, corresponding to subgraphs with 
g = 2, 3, 4, 5 and 6 edges. In each family, the scaling is 
the same, but the prefactors generally differ and depend 
on the dimensionality and on the form of F(r). 



D. Geometric networks with arbitrary degree 
sequences 

Real- world networks often have degree sequences which 
are quite different from Poissonian. For example, many 
networks have heavy-tailed deg ree distributions [ESS 

EL El El El IH M El EI . 

A heavy-tailed degree distribution in random networks 
has been shown to strongly affect the counts of certain 
subgraphs ,44]. In particular, certain subgraphs appear 
much more often in random networks with heavy-tailed 
degrees than in Erdos networks. For example, subgraphs 
in networks with a scale-free out-degree P(k) ~ fc -7 and 
compact in-degree have been shown to scale with network 
size N as 0: 



(G) ~ A^ Q 



(13) 



where a is related to the number of subgraph nodes n, 
subgraph edges g and maximal subgraph out-degree s: 



n — g + s — 1 
a = {n — g + s — 7+I 
n-g 



7 < 2 
2 < 7 < s + 1 
7 > s + 1 



(14) 



The geometric models discussed above have a Poisso- 
nian degree distribution. An interesting extension of ge- 
ometric models, which allows for heavy-tailed sequence 
has been recently studied by Havlin and colleagues. The 
scaling of path lengths in that model was derived 39] . 



Here, we consider subgraphs in a related, directed lat- 
tice model, with an arbitrary outgoing degree sequence 
P(k). In this model, each node in the lattice is assigned 
an degree k drawn from the distribution P{k). k out- 
going edges are then randomly connected to other nodes 
according to the connectivity function F(r). This results 
in a geometric model, with outgoing degree distribution 
P(k), which can be heavy-tailed, and a compact incoming 
degree distribution. We note that several real-world di- 
rected networks have compact in-degree and heavy tailed 
out-degree, including biological [3 US an( A technological 
networks |54ll55|. 

We now derive scaling relations for the number of sub- 
graphs in geometric networks with a heavy-tailed outgo- 
ing degree sequence: 



P{k) ~ k^ 1 . 



(15) 



We consider the limit where the hubs in the network do 
not exceed the typical size of the neighborhood of each 
node, which scales as R d . The mean hub size scales as [H 



(16) 



Thus, we consider the case where N 1 /^^ « J? d « JV. 

The network can be considered as a collection of N/R d 
subnetworks of size N n = R d , each with scaling according 
to Eq. 13. For the entire network, one finds that the 
number of subgraphs G scales as 



(G) - (N/R d )N n a = NR^- 



(17) 



with a given in Eq. 14. All subgraphs numbers scale as 
N, and have an R dependance that depends on 7 and 
the subgraph topology. Subgraphs with large maximal 
out-degree s tend to appear more often than subgraphs 
with smaller s. 
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field strength (1-p) field strength (1-p) field strength (1-p) 

(a) ~ (b) (c) 



FIG. 4: Three subgraph ratios for different field strengths in an anisotropic geometric model with gaussian connectivity function. 
d = 1,N = 2000, R — 30, (k) = 8. Shown are simulation results (a-c) and theoretical curves, obtained from Table II (a) and 
Appendix B (b,c). Numerical standard errors are smaller than the dot size. 



E. Invariant ratios 

It is of interest to find invariant measures which apply 
to geometric networks independent of the dimensionality 
d and the form of F(r). Several such invariants can be 
found. An example is the ratio of the numbers of two 
subgraphs, the feedback loop and the feed-forward loop 
(subgraphs 6 and 7 in Table II). The ratio of these two 
subgraphs is 1 : 3, regardless of d and F(r). Similarly, 
the ratio of subgraphs 9,10 and 11 is 1 : 2 : 1. The 
ratios of subgraphs 1,2 and 3 is 1 : 2 : 1. Table III 
shows these ratios (ratios 1 and 2) for several real world 
networks. The ratios in real world networks generally 
differ significantly from the ratios in isotropic geometric 
networks. 

The above ratios are invariant in the poissonian geo- 
metric model, but can change if one considers geometric 
networks with an arbitrary degree sequence. For exam- 
ple, the feed-forward loop subgraph (subgraph 7 in Table 
II) has a node with 2 outgoing edges (s = 2), whereas the 
feedback loop (subgraph 6 in Table II) has s = 1. This 
would result in an increased abundance of feed-forward 
loops for geometric networks with heavy-tailed degree se- 
quence. 

We therefore sought ratios of subgraphs which do not 
depend on the degree sequence. These are subgraphs 
which have the same subgraph degree sequences. The 
class of four-node directed subgraphs contains several 
examples of subgraphs with identical degree sequences. 
Two examples are shown in Table III (ratios 3 and 4). 
These ratios are invariant with respect to the degree se- 
quence of the network model, as well as to the dimen- 
sionality d and the form of F(r) 



F. Directed anisotropic geometric model 

In some systems there are preferred directions in which 
the probability of a directed connection is larger than in 
other directions. To model this, we consider a d = 1 



lattice, and introduce a bias that favors edges in one di- 
rection. In this probability of a directed edge to 
the right is F ri gi lt — pF(r) and an edge to the left is 
Fi e ft — qF(r), such that p + q = 2 (the sum=2 is cho- 
sen to preserve the same (k) as for the isotropic model). 
The calculation is somewhat more intricate in this case, 
as the overlap integrals must be evaluated for different 
orderings of the nodes (Appendix B). The expression for 
each subgraph is multiplied by a 'field factor', which de- 
pends on the subgraph topology. The results for 3-node 
directed subgraphs are shown in Table II. 

The scaling is the same as in the isotropic case, except 
in the limit p = (or q = 0). In the latter limit, no cycles 
are allowed (the network is a directed acyclic graph) . The 
only closed three-node pattern is the feed-forward loop 
(subgraph 7 in Table II). 

The relative abundance of the various subgraphs de- 
pends on the 'field' p (Table II). Subgraphs with mutual 
edges (edges in both direction between two nodes) are bi- 
ased against since mutual edges always contain one edge 
that goes against the field. The ratios that were invariant 
in the isotropic case are no longer independent of d and 
F(r). Several of these ratios are plotted as a function of p 
in Fig 4. Table III shows the ratios in real-world networks 
for the fully biased anisotropic geometric network. 



III. DISCUSSION 

This study presented analytical results for the sub- 
graph content of geometric network models. We find 
scaling rules for the number of appearances of each sub- 
graph as a function of network size and lattice dimension- 
ality. The scaling is very different from Erdos-Renyi net- 
works for most subgraphs. We find certain ratios of sub- 
graph appearances in isotropic directed geometric net- 
works, which are 'invariants' in the sense that they do 
not depend on dimension and connectivity function. 

Geometric networks show distinct network motifs. All 
of the subgraphs scale as N 1 , whereas they scale as N n ~ 9 
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in the corresponding random networks with the same de- 
gree sequence. Therefore all subgraphs with g > n are 
network motifs in geometric networks. In most real- world 
networks studied so far, only a, subset of these subgraphs 
are network motifs @, Q, 0> E3 > suggesting that ad- 
ditional constraints and optimization is at play, beyond 
isotropic geometric constraints. 

The network motifs in real-world networks appear to 
be 'extensive variables', in the sense that their concen- 
tration c = G/N does not decrease with N, but scales 
as N° 25j, whereas their concentration in correspond- 
ing randomized networks decreases with N. A similar 
property is found for the concentration of motifs in ge- 
ometric networks, which scales as c = G/N ~ N°, and 
c = G/N ~ TV"" 5-1 < TV -1 in randomized networks. 

Table III shows that the abundance of feedback loops 
relative to feed-forward loops in social and world-wide- 
web networks is much less than expected from an 
isotropic geometric constraint |6fi| . 

Social networks which represent positive sentiments 
between individuals in a gro up are known to be rich in 
transitive relations 8, 59, 60]: If X 'likes' Y and Y likes 
Z, X tends to also like Z . In these networks, intransi- 
tive triplets of nodes (X — > Y, Y — > Z, but X Z) are 
known to be relatively rare. Feedback loops might be rare 
because they contain 3 intransitive triplets and no tran- 
sitive triplets. The feed-forward loop, on the other hand, 
contains 1 transitive triplet and no intransitive triplets. 
This might also explain the relative rareness of subgraph 
10 of Table II, which contain 1 transitive and 2 intransi- 
tive triplets, as opposed to 2 transitive and intransitive 
triplets in subgraphs 9 and 11. Similar transitivity may 
characterize WWW links Q . 

In the neuronal network of C. elegans feedback loops 
and feedback loops with one mutual edge (subgraphs 6 
and 10 in Table II) are much less abundant then expected 
in isotropic geometric networks (Table III). We analyzed 
two versions of the neuronal network - a network which 
includes all synaptic connections, and a network which 
includes only 'strong' connections, where neurons are 
linked by an edge only if they have 5 or more synapses 
connecting them. 

Interestingly, when one considers subgraphs with no 
mutual edges (ratios 1, 3, 4 in Table III), the full neuronal 
network of C. elegans displays similar subgraph ratios to 
a highly anisotropic geometric model (with field p<l). 
Such a field may represent a directed flow of informa- 
tion from sensory neurons, generally located in the head, 
to motor neurons. This model cannot however explain 
the relative abundance of mutual edges, which would not 
appear in highly biased anisotropic geometric networks. 

The subgraph ratios in the neuronal network of strong 
connections is not consistent with a highly anisotropic ge- 
ometric model (Table III) . It seems that the motifs in this 
network do not solely stem from geometric constraints, 
and that additional optimization constraints based on bi- 
ological functionality are at play |sl ITTl HI 111 HI |25| . 
Indeed, in transcriptional biological networks, network 



motifs have been experimentally shown to function as 
pulse generators |2]j , asymmetric filters j2^ , response ac- 
celerators [2^ and temporal pattern generators [3(], . 

The present scaling results also apply to the class of 
small- world network models Q. These models are ob- 
tained by rewiring a limited fraction of the edges which 
leads to a substantial drop in the network diameter but 
in a small change in the clustering coefficient. 

The geometric model can be extended in many ways. 
One can assign an arbitrary clustering coefficient for each 
node, in addition to an arbitrary degree, by specifying 
both the number of edges ki and a neighborhood Ri for 
each node (Eq. 10). This can form a more stringent 
comparison to networks with broad deg ree sequences and 
clustering coefficient sequences 0, . In addition, bi- 
ases towards mutual edges can be added, in addition to 
the anisotropy field discussed above. 

It would be of interest to extend this study to under- 
stand the general relationship between global constraints 
on a network and its local structure |6lL 1621 l63l |6J| . 
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APPENDIX A: OVERLAP INTEGRALS FOR 
THE HARD-CUBE CONNECTIVITY FUNCTION 



-R 



y R 



FIG. 5: Overlap area 0(y) = 2R — y shaded in gray. All 
lattice points within this area can connect to both the origin 
(0) and to y. 



The intersection integrals reduce to a simple form in 
the case of the hard-cube connectivity function F c . As 
an example we will consider triangle subgraphs: 



(Ga) 



N 



N 

7T 



F{0,y),F(0,z),F(y,z)dydz 



(Al) 

Where Id is the d-dimensional overlap integral. In the 
hard-cube case it is sufficient to calculate the integral 
in one dimension I\. For this connectivity function the 
axes are separable, and the extension to d dimension is 
straightforward [53, |58j : 



Id = h L 



(A2) 
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I\ reduces to a one-variable integral 



( k 



h= / F(0,y),F(Q,z),F(y,z)dydz-- 

0{y)dy=> 



( k ^ 



\(2R)JJ_ R vy/ » \(2R)y. 

/ k 



f (2R-\y\)dy = 

J-R 



\(2R)J ( 3 ^ = ^ 1 ( A3 ) 

Where 0{y) is the overlap area, which includes all posi- 
tions of z which are in the range R of both and y (Fig 
5). 

Extending these results to all directed types of trian- 
gles (Table II) and arbitrary dimension yields: 



(G) = a ( — ) NR^ n - g - 1)d k 9 



(A4) 



where a is a symmetry function equal to one over the 
number of permutations of nodes that lead to an isomor- 
phic subgraph. Eq. A4 contains the scaling of Eq. 9. 
The prefactors depend on the symmetry function and on 
the dimension and decrease with d. The extension of Eq. 
A4 to larger subgraphs is : 



(G) = a (^J NR^-a-^k 9 



(A5) 



where vq is the factor stemming from the Id overlap 
integral, which depends on the non-directed version of G. 
These factors have been calculated in [5]J for subgraphs 
of up to 7 nodes. The prefactors for all non-directed 
4-node subgraphs is given in Table II. 



differently. An example for the feedback loop is shown in 
Fig. 6. For this subgraph symmetries lead to a simplified 
expression for the field factor ip: 



p = 7;pq(p + q) =pq 



(Bl) 



For n-node subgraphs one should calculate the n — 1 di- 
mensional integral over n\ hyper- volumes, corresponding 
to all linear ordering of the subgraph nodes. For the case 



q 2 p 

y z 


PQ / 

y z X 

/ q 2 p 

*- ► 


/ 

X z y 


p 2 q 

z y 



FIG. 6: Orientations used in the calculations of the field factor 
for the feedback loop in the anisotropic geometric model. The 
y — z plain is divided into 6 regions of integration, and the 
weighting factors are shown above the subgraph diagrams. 



APPENDIX B: 



CALCULATION OF FIELD 
FACTORS 



Here we describe the method for calculating the field 
factors for the three-node subgraphs in the Id anisotropic 
geometric model. We fix one node at the origin (0). The 
two remaining nodes are denoted by y, z. The formulas 
for the subgraph occurrences are similar to Eq. 2, but 
here different integration domains are weighted by dif- 
ferent field factors, depending on the orientation of the 
three nodes. The y — z plane is divided into 6 regions 
which correspond to the 3! linear ordering of the three 
subgraph nodes. These regions are generally weighted 



of the 4-node subgraphs of the invariant ratios presented 
in Fig. 4(b,c), this calculation results in the following 
theoretical curves: 

b(p 4 + q 4 ) + 4(a + c)pq( P 2 + q 2 ) + 6bp 2 q 2 
n a(p A + q 4 ) + 2bpq(p 2 + q 2 ) + 2(a + 2c)p 2 q 2 1 ' 

where a = 1.1 • 1(T 2 , b = 8.8 ■ 1CT 3 , c = 5.4 • 1(T 3 , and 

1 h(p 2 + q 2 ) + (2e + 2f + g- h)pq 



2 e{p 2 + q 2 ) + (2/ + g + h)pq 



(B3) 



where e = 2.35 • 1CT 3 , / = 1.9 • 1(T 3 , g = 3.8 • 1CT 3 , h = 
1.8 • 1CT 3 . 
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